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1 MTHRTL.SRC JMTHGEXP.MAR; 1 
-TITLE MTHSGEXP 5 geting Exponential Function 
IDENT /1-006/ ; : MTHGEXP.MAR Edit :RNH1006 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
eoRPORAT on NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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fACILITY: MATH LIBRARY 
ABSTRACT: 
MTHSGEXP is a function which returns the G floating point 


exponential of its G floating point argument. 
The call is standard call-by-reference. 
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VERSION: 1 
HISTORY: 
AUTHOR: 
Steven B. Lionel, 15-Jan-79: Version 1 
MODIFIED BY: 
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-SBTTL HISTORY ; Detailed Current Edit History 
Edit History for Version 1 of MTHSGEXP 
BH = Adapted from MTHSGEXP version 1-008. SBL 15-Jan-79 
17004 


MTHSGEXP . Floatin 
1 fro Y; De 


¢ - Corrected a typo in the e. JBS ; -JUL-1979 
Use only through R6. SB 1SSeptei9 
passe large argument logic to ners lose ay agit asence in 


EMOD - # ogo greater tha Te 
SSIGNAL al MTHSSJACKET_TST 


-005 - RAN GSesep on calls. to MTH 
-006 - Elim conen ecfoBs short Literals 


MTH 
1 


000407¢ 
0001760 


hoya yon + 4 


SVU" WOOO 


nm 
oO 


OVOMoO CFP OMOOOO- 


Nomen RS 
POW 


i +0 5 

oro» 
wormor 

eaomons 


esouSS oe HOLL LOL OLWEOLOS 


ony 
om 


BE 
sm 


Ww 


Dan SAAS tote Bee he ee ete ee eee Ee NS este. 


DWONANE WIN @§ O OD NOUN EWN 9 CONAN 3 OOD NOAU SWAN (OOO NOU WA) OUWOONOu 


BB a td te 8 ed at tS et 
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$ INCLUDE FILES: 
EXTERNAL SYMBOLS: 


é .OSABL 
"EXTRN 


~EXTRN 


-EXTRN 
oEXTRN 


ACMASK 


X_16L0G2E = 


MACROS: 
SSFDEF 


»PSECT 


: OWN STORAGE: 
: CONSTANTS: 
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EPEC or noaw RESTON SEEDY ARATE SEGRRB eens Mo 


DECLARATIONS ; Declarative Part of Module 
MTHJACKET.MAR 

GBL :; Declare all EXTRNs explicitly 

MTHSSSIGNAL 3 SIGNAL “ven error 

MTHSSJACKET_TST 3 vest to see if called with CALL or 

MTHSK_FLOUNDMAT ; Underflow error code 

MTHSK_FLOOVEMAT ; Overflow error code 


EQUATED SYMBOLS: 


= “M<IV, Re R3, R4, RS, R6> ; ENTRY mask + int ovf enable 
“001354 6 


546 ; Extension for operand in EMOD 


; define SF$ (stack frame) symbols 


PSECT DECLARATIONS: 


_MTHSCODE PIC, SHR,LONG, EXE ,NOWRT 
3 program section for math routines 


none 


; Table to be used for scot ing. These constants here have been 
er 


; verified by M. Payne using 
*0040020,0 
0,0 


TABHI: .WORD 
ow RD 


program ROOT16 on PDP-10. 
2**(0/16) = 1.0 
3 2e0(1/16) 
2**(2/16) 
3 2e0(3/16) 
3 20#(4/16) 


239 
i 
i 3 2e2(5/16) 
é 3 2%*(6/16) 
iH 
i 


3 2e2(7/16) 
3; 2*(8/16) 
3 202(9/16) 
2**(10/16) 
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WORD *0040031,°0142221 : 20#(11/16) 
; woe comma stage): eect 
“WORD — *0040032,°016 : 2#8(12/16) 
WORD %0114330,°0151685 
WORD 0040054, 20016 : 2#(13/16) 
. “0026035. -0084030 2*#(14/16) 
7 . 20136 **0 3 . 
? "WORD *0040036,40129957 : 20*(15/16) 
? : 90121264,°0110331 


\\N NOTE!!!: The decimal equivalents in TABLO are taken 
directly from MTHSDEXP, They do not correspond to the 

actual value of the G floating equivalent. When an entry 

in TABHI and its pal ag oti entry in TABLO are added, the 
should equal the correct fractional power of 2 to 74 bits. \\\ 


TABLO: eee 8 3 DECIMAL: 0.00 
“WORD *6036270,*0123056 ; DECIMAL: 0.2252169616881804D-17 
“WORD *6036317,°0115760 ; DECIMAL:-0.2712242510500122D-17 


“WORD 
34 "WORD 

09 = WORD “6036271,°0130177 : DECIMAL: 0.5861402647731367D-17 

OA ORD 30036246,*0172153 3 DECIMAL: 0.1206457647223494D-16 

OAS -WORD *6036252,°0155011 3; DECIMAL :-0.8930877995013540D-17 

oad -WORD *6036235,“0061630 ; DECIMAL :-0.2373071989573779D-17 
a : WORD 85036266,*0031114 ; DECIMAL :-0.6257240830881880D-17 

06 “WORD 85036306,*0010546 3: DECIMAL :-0.1340620676392399D-16 


RD 0 

“WORD *6036316,°0175242 > DECIMAL :-0.7084371812598154D-17 
-WORD *6036266,°0164761 : DECIMAL :-0.3768379065187162D-17 
*WORD 85036234, *0076107 : DECIMAL: -0. 304838430961 3603D-17 
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D8 
43 -WORD 0 
2 — 9° 36267 ,*0120715 : DECIMAL:-0.12766242353000400-17 
4 -WORD 0 36241 ,*0010146 ; DECIMAL: 0.1845830375854930D-17 
i -WORD “9 36242 ,°0166403 3; DECIMAL: 0.5075495866202897D-17 
F8 ; WORD *6036300,0130756 : DECIMAL: 0.48228430606756190-17 
F -WORD 0, 
: Constants used in evaluation of polynomials - small arguments 

GXPTB1: .WORD *0037432,°0001131 
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xponential Functi 


tion y 
; Declarative Part of Modu 6- 
WORD 


GXPTAB: rm 


“WORD 0 
GXPLN=<.-GXPTAB>/8 


6_16L0G62_E: 
“WO 
“WO 
G_LN2_0V_16 
7QU 


G_LN2_0V_16 
QU 


HI: 


LO: 


ah 
*0 
*0 
*0 
*0 
*0 
*0 
*0 
*0 
*0 
*0 
4 
4 
7 
4 


~ 0 

0 

™ "0177777 
*01 "0177777 
PA 0 

o8 


; Constants used in evaluation of 


pines cot 
Sa 
SHES) 
70037065.0013617 
sie 1220124727 
*0037 33° 0131253 
STE: att 
eT 
“017780 2014261 
FE RS 


“0062455,°0101376 


“XOOO0OF EFA2E423F C6 


“X3B3ABC9EF 79A3D5C 
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DECIMAL: 


; DECIMAL: 


DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 


DECIMAL 


0.2480427857745020D-04 
0.1984369200268758D-03 
0. 1388888879690042D-02 
0.8333333262370290D-02 
0.41666666666679500-01 
0. 1666666666667437D-00 
0.5000000000000000D-00 
0. 10000000000000000+01 
0. 10000000000000000+01 


no. of entries in table 


polynomial - regular args 


DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 


DECI 


0.30771307094302400-15 
0.5682419384166091D-13 
0.9181219559808114D-11 
0.1271587192556359D-08 
0.1467610032291993D-06 
0.1355080777949815D-04 
0.93838479280898720-03 
0.4332169878499658D-01 


MAL: 0 
no. of entries in table 


; LOG2(E) * 16 


Hi 39 bits of ln2/16 


; Low bits of ln2/16 


Page (3) 


to the signal mechanism vector CHFS$L_MCH_RO/R1 by LIBSSIGNAL). Associated 
message is: "FLOATING OVERFLOW IN MATA LIBRARY’. Result is reserved operand 
-0.0 unless a user supe (Ned (or any) error handler changes CHFSL_MCH_RO/R1. 

MTHS_FLOUNDMAT if X =< -709 and caller has hardware enable set. The 
repute et set to +0.0. Associated message is: ‘FLOATING UNDERFLOW IN MATH 


NOTE: This procedure disables floating point underflow, enable integer 
overflow, causes no floating overflow or other arithmetic traps, an 
preserves enables across the call. 
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MTHSGEXP - Sranderd G Floating EXP pas et 94335 :99 MTHRTL.SRC JMTHGEXP.MAR; 1 = (4) 
. ; «SBTTL MTHSGEXP = Standard G Floating EXP 
A j p+4 
A 3; FUNCTIONAL DESCRIPTION: 
“ 33 : EXP -G floating point function 
i ; Uses a Chebyshev approximation, with overhang on Last step. 
A 43; 
. rf: : CALLING SEQUENCE: 
n ¢$ ; Exponential.wg.v = MTHSGEXP(x.rg.r) 
A 48 : INPUT PARAMETERS: 
00000004 A 50 LONG = 4 ; define Longuord multiplier 
00000004 ~ 2) x = 1 * LONG ; Contents of x is the argument 
r 38 IMPLICIT INPUTS: none 
: 33 OUTPUT PARAMETERS: 
As 33 VALUE: G floating exponential of the argument 
As 23 IMPLICIT OUTPUTS: none 
. ) SIDE EFFECTS: 
ms e§ Signals: MTHS_FLOOVEMAT if X > 709 with reserved operand in RO/R1 (copied 
2 
7 
3 
5 
4 
78 
74 


7 

7 

78 ENTRY MTHSGEXP, ACMASK 3 standard call-by-reference entry 

79 3; disable DV (and FU), enable IV 
MTHSFLAG_JACKET ; flag that this is a jacket procedure 


6D 00000000'GF 9% MOVAB G*MTHSSJACKET_HND, (FP) 


3; set handler address to jacket 
handler 


in ‘ se of an error in special JSB 
: ne 

; RO/R1 = user's ore 

; RO/R1 = special EXP(RO/R1) 
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3; return - result in RO/R1 
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5". Special GEXP routine aise 94:55:99 EATHRTL. SREINTMGEXP MAR; 1 _ ) Ta 
-SBTTL MTHSGEXP_R6 = Special GEXP routine 


Special GEXP = used by the standard, and direct interfaces. 


CALLING SEQUENCE: 
pove anything needed in RO:R6 
JSB MTHSGEXP_R6 
return with result in RO/R1 


Note: This routine is written to avoid causing any integer overflows, 
floating overflows, or floating underflows or divide by 0 conditions, 
whether enabled or not. 


REGISTERS USED: 
sha - floating argument, then result 
RS = temp 


oat 


- GFL 
MTHSGEX 


MTHSGEXP 
1-006 


ve 
xz. 


: input in RO 


Re/ 
RS - integer scratch 
R6 = integer part of X * LG2(E) * 16 (16I+4J) 


MTHSGEXP_R6:: ; special GEXP routine 
MTHSGEXP_R?:: 3; Release 1 name 


The preliminary test for overflow works as follows: First, the sign bit is 
cleared leaving the first word of the ‘Xi. Then, 1024-4 (bias-4) is sub- 
tracted, leaving an exponent biased by 4 in bits 14:4 and the first four 
fraction bits in 35:0. This rebiased value is compared against 230 (decimal). 
The comparison can have 3 outcomes. If the rebiased value is now negative, 
this means that the true exponent is < -4 - gris is a BLSSU test. If the 
rebiased value is posts ive but greater than 230 (decimal), then the ‘Xi is 
greater than or equal to $6, which is pyerentece overflow or underflow, 
epending on the sign of X - this is a BLSS test. Otherwise, X is somewhere 
in the range for the standard evaluation, and flow continues. 


Sete Ge Ge Ge Se Ge Se Se Se Se Gee 


52. 50 8000 8F 


53 26 3FCO 8F = AB 
5 QOE6 8F_ BI 
os 


BICW3 #*x8000, RO, R2 
SUBW3 «= #*X3FCO, R2, R3 
CMPW 0 #*XE6, RB 
BLSSU  SMTST 


Preliminary test for over/under flow 
R2 = exponent bits only 

R35 = 4 + unbiased exponent 

unsigned compare of :Xi with 732 

to more tests if LSSU 

else, -4 < unbiased exp < 11 

no exceptions in EMODG or APPROX 


Check for loss _of_significance in 
EMOD ( IX! >= 2887 


4080 8F 52 B61 CMPW R2, #*xX4080 


So 
=] 
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HLH FAA NIAINIMIAIGEIPIPINPINPIPUDINPININID) 2 I | MOQ OODOOOOOCOOOOOO0O0O0O0O0OMoM°§9 
WRI § OS OD NAW EWN 9 OO NAME WIN © OD NAU EWN CE ODNAU EW O OONOAUES WOOO 


28 «19 BLSS EVAL ; No loss of significance 
SIX} >= 2004 FRQO will lose significance so the interger and fractional 
; parts of X*16/[n2 must be obtained in seperate steps. 
52. 50 89 AF 45FD : MULG3 6_16L0G2_E, RO, R2 ; Get integer part of X*16/Ln2 in 
28 2¢ 4AFD CVIGL R2, R 3 R6 (=I+J) as a longword and in 
2 6 ra CVTL R6, R ; R2/R3 in G rormat 
54 52 B3 AF 45FD MULG G_LN2 OV_16_HI, R2, R& ; Get frac Ma part of X*16/\n2 = 
50 54 42FD SUBG Ra, RO 3 16/Lin2eC KX = (14J)*bn2/16 J 


D 12 
ponential Function 


MTHSGEXP G Floating Ex 16-SEP-1984 01: AX/VMS M v04-00 P 
1-006 ATHSGEXP.R6o= Special GEXP routine ore 1 obe te $f: 7 YRTART eae SaTH HGEXP.MAR; 1 oo (3) 
52 so” af c3tD iF te SUBS G_LN2_OV_16_LO, R2 ; in RO/R1. 
50 98 AF 44D O1F 4 MULG2 G_{ R : 
OA 11 i ‘3 BRE ty ey § 
50 56 50 1760 8F 91 AF S4FD OF i EVAL: EMODG G_16L0G2_E, #X_16L0G2E, RO, Mo Me LEaCe) seh 
; ge w 
o i 3 ante er part in R6 (=161+J) 
§ : fraction in RO/R1 
FF3C CF O08 50 SS5FD 27 APPROX: POLYG RO,#GXPLN-1,GXPTAB ; use Chebyshev series 
C 22 3; with last coefficient 0 
bc 28 3 $o that Last ADDG has overhang 
55 56 FFFFFFFO 8F cB 0c 58 BICL #-16 Ro RS 3 RS 2 J 
50 EDEG CF45 44FD 0214 359 MULG2 TABHICRSJ, RO : Ps MUL oP 2e #64716) 
50 FESF CF45 4OFD 1 60 ADDG TABLOLR5J, RO 3; add in of 2**(J/16) 
50 FDD8 CF45 4OFD 8 2 ADDG TABHICRSJ, RO 5 ond Ey HI of 2%*(J/16) 
56 OF CA 0229 86 BICL #15, R6 3 = b 
08 13 022C 64 EQL 0$ : Re 1=0, then done 
50 56 CO O22E 65 ADDL2 Rk6, RO 3 Add I i exponent. 
0231 06 3 MUL by 2**1 by exponent addition 
OF 50 61 0231 6 CMPW a*XF 3; test for over/under flow 
40 15 8 4 68 BLEQ RO cERT 3 see what exception is if neg or = 0 
05 $ $ 208: RSB 3; otherwise return result in 
1% 19 : ‘teeoagt ea tion if exp+4 > 14 
3: exception exp+4 > 
3ca0 ae 83st 0930308 CNPW RZ, #*X3C80 ; eliminate underflow. from APPROX1 
08 19 O23E 37% BLSS 108 : bypass if E**ARG = 1 
0540 376 s+ 
8 40 5S ; Use Chebyshev series for small arg 
$500,379 
FEB9 CF 08 50 55FD QO rb 9 POLYG RO,#GXPLN1-1,GXPTB1 Use Chebyshev series 
024 1 : last term is 1; this will 
47 6 : give desired overhang. 
05 rt 4 ? RSB 3 answer is OK, return 
50 08 SOFD iB 5 10$: MOVG $*#1, RO 3; Ex*X is 1, store it 
05 4C § RSB 3; and return 
40 
Be 
rt y ; Handlers for software detected over/underflow conditions follow 
50 53FD 0240 392 20$:  TS1G RO : if big ARG > 0 goto OVERFLOW 
38 AR O58 $8 BGEQ OVER . . 
32 ; Underflow; if user has FU set, signal error. Always return 0.0 
39 GNDER 
5 dc 4 MOVPSL R : Re = us er’ or jeghet routine's PSL 
00000000 ' GF 2 4 4 9 CALLS #0, GRTHOSJACKET 1ST 3; RO = ag « from jacket routine 
04 E 5B 400 BLBC RO, 108 > branch ar did JSB 
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1-006 MTHSGE ENP Ree Spec al GEXP routine g- SEP=1984 t 9] 3h: % YATHRTL. SRESATHGEXP MARS 1 . (8) 
52 O4A C 5 4 MOVZWL SFSW_SAVE_PSW(FP), R2 ; get user PSL saved 
8 3 és 4 10S: CLRO RO ~ ni : i = result 5 /AtSESTCMAL wi will save in 
4 3 CHFSL_MCH_ Ro /R1 so ny handler cen ee 
0D 52 06 €1 64 404 BC R2, 208 3; has user enabled tloa ing underfl 
7 bp 68 405 PUSHL (88, i yes, return PC from s spec al rout in 
7E 00°8F A A 4 § MOVZBL #MTHSK_FLOUNDMAT, -(SP) ; trap code for hardware f underflow 
é ¢ 3 convert ¢ to. pb FLOUNSMAT’ H jn it VAX-11 
3 exception 
00000000'GF 02 Fe 3 4 § CALLS #2, G*MTHSSSIGNAL 3: si nal tins PC) 
0 ay 20$: RSB : return 
8 tig EXCEPT 
56 6S 7 4) TSTL R6 ; test sign of I; if 1 <0 
dB «19 4 aie. BLSS UNDER 3; go to underflow handler 
gh $18 : 3; Signal floating overflow, return reserved operand, -0.0 
7A 418 dver 3; else process for overflow 
g OD 7A 419 PUSHL (SP) 3; return PC from special routine 
7E OO°8F QA 7C «420 rt #MTHSK_FLOOVEMAT, -(SP) ; hardware floating overflow 
5 1 #O 79 0280 421 ASHQ #15, #T, RO : RO/R1 = reeus = reserved. eperend 
0 He 4 § ; “0.0. RO/O will be copied t 
8S 4 4& 3 signal mechanism vector (CHESL -MCH_RO/R1) 
0284 424 3 so can be fixed up «i any error 
0284 425 ; handler 
00000000'GF 02 Fe 0 ee 426 CALLS #2, G*MTHSSSIGNAL 3 si + i (condition, P 
05 0288 427 RSB 3 return = RO phe on By from CHF$L_MCH_RO/R1 
028C 428 
bsec $93 
028C 430 .END 


ee ne ee 


F 
MTHSGEXP 3; G Floating Exponential Function 


Symbol table 


ACMASK = 4, 
APPROX 


DwDwD 
rorore 


oa 
= 
v 
z 
_ 
“un 
PoOoorooyn"on~ 


22 
— 
nee 

m 

< 

v 

a 

» 
oo 
ooo 
ooo 
ooo 
ooo 
tb 
ow» 
~o 
2wwD 
aa 

aed 
SOOO SCOOOOOOSOO COCOO ooo 
MMMM WOOMRMmMOCr 


1 

0 

ND eteereee x 

ST steenere x 
geeteene 


eeeeeere x 
keeeeere x 


WT HSK. Ff LOUNDMAT 


ov 0000027A R 
Oy eu SAVE_PSW = 00000004 
SMTST si 7R 
TABHI 00000 R 
TABLO 0000 3 R 
UNDER 0000 R 
x = 00000004 
X_16LOG2E = 00001760 
Serr ana enon sasome + 
! Psect synopsis ; 
PSECT name Allocation PSECT No. Attributes 
ABS . 00000000 < 0.) 00 ¢ QO.) 
S$ 00000000 ( 0.) 4 ( 1.) NOPIC US 
MTHSCODE 0000028c ¢« 652.) 02 ¢ 2.) PIC US 
(uveneninenninnmmamnmnel 
! ; Performance indicators ! 
Phase Page faults CPU Time aateses Time 
Initialization 3} :00: 8-1) 700. jf 
oe processing 198 :00: 9-88 Ree 
pyane’ table sort 99 : 0:03 
od, table output 4 ° 00.04 : 
Psect synopsis output :00:00.01 ‘ 
Cross-reference output 3: :00. ! : 
Assembler run totals 364 90:09 


[oe working sgt Limit was 900 page 
2 bytes (17 pages) of virtest aonere were used to buffer the intermedi 


NOPIC USR CON ABS 


"S-SEp=1986 91:55:97 EATMNTE.SREaRTMGExPcman:1 | "#9 


zz 
o 


LCL At he NOEXE NORD ey NOVEC pa 
R CON ABS LCL NOSHR EXE RD T NOVEC BYTE 
R CON REL LCL SHR EXE NOURT NOVEC LONG 


ate code. 


11 
(5) 
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MTHSGE XP : G Floating Exponential fmatton ” i a + P= 13 91: $f: 99 ¥ 


VAX=-11 Macro Run Statistics P= 


bs re were 9 pages of symbol — space ot d to hold 56 non-local and 5 local symbols. 
source Lines were read in Pass 1, produc +o ob ect records in Pass 2. 
3 pages of virtual memory were aed to Poet ine macros. 


Proce nmemooas er ascamoen sass wmanwa} 


! Macro Library statistics ! 


$eoreaeecer one tense n creme mocen $ 


Macro Library name Macros defined 
“$255$DUA28: CSYSLIBISTARLET. MLB; 2 4 
88 GETS were required to define 4 macros. 


There were no errors, warnings or information messages. 


AX/VMS Macro V04-00 
MTHRTL.SRCIJMTHGEXP.MAR; 1 


ro Ht 


MACRO/ENABLE=SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHGEXP/OBJ=OBJ$:MTHGEXP MSRCS:MTHJACKET/UPDATE=(ENHS :MTHJACKET) +MSRCS: 


NT CORPORATION 
D PROPRIETARY 


AH-BT13A-SE 
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